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I We study radiation-matter interaction in a system of ultracold atoms trapped in an 

optical lattice in a Mott insulator phase. We develop a fully general quantum model, and 
we perform calculations for a one- dimensional geometry at normal incidence. Both two- 
and three-level A atomic configurations are studied. The polariton dispersion and the 
reflectivity spectra are characterized in the different regimes, for both semi-infinite and 
_ finite-size geometries. We apply this model to propose a photon energy lifter experiment: 
O ! a device which is able to shift the carrier frequency of a slowly travelling wavepacket 
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without affecting the pulse shape nor its coherence. 
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^ '. 1 Introduction 

■ It is a general fact of electrodynamics in continuous media that matter acts as a sort 
of potential for light: the dielectric constant in Maxwell equations plays in fact a role 
mathematically equivalent to the external potential in the Schrodinger equation [1]. In 

O ! particular, remarkable features appear for light propagation through a system with a 
periodic modulation of the dielectric properties on the scale of the radiation wavelength 
^ ' [2,3]: in analogy with electron propagation in crystalline solids, these structures are then 

■ called photonic crystals (PCs). According to the Floquet-Bloch theorem, the discrete 
translational symmetry garantees in fact the conservation of the Bloch wave vector in the 
first Brillouin zone, the spatial periodicity of the eigenf unctions in the elementary cell, 
and the separation of the energy spectrum in bands and forbidden gaps [4,5]. Depending 
on the frequency dependence of the dieletric response of the material, PCs are generally 
classified into two broad categories: passive, if the underlying media are dispersionless, 
and resonant [6, 7]. 

An interesting example of resonant PC can be realized using ultracold alkali atoms trapped 
in the periodic potential of an optical lattice [8,9]. In particular, a mile-stone in this field 
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has consisted of the achievement of a Mott Insulator (MI) phase [10, 11], in which an 
integer and constant number of atoms are trapped at each site of the periodic optical 
potential formed by the interference pattern of several laser beams. If temperature is 
low enough, atoms are frozen in the lowest vibrational level of each well and the system 
periodicity is extremely regular: almost no impurities are in fact present (i.e. missing or 
extra atoms), and there are no phonons. Thanks to the simple level structure of alkali 
atoms, one can selectively address specific transitions, so to realize e.g. two or three- 
level models. From this point of view, the MI can be seen as an extremely resonant PC, 
where the Bragg scattering processes due to the periodic arrangement of atoms have a 
strong interplay with the atomic optical resonances. The cleanness of the system, and 
the weakness of spurious non-radiative effects guarantees that coherence can be preserved 
for very long times during radiation-matter interaction. This is crucial to manipulate the 
photon propagation without losing its coherence. 

In this paper, we develop a general and fully quantum theory to describe the radiation- 
matter interaction in these systems. Although the model is perfectly general, we specialize 
our analysis to the ID case with normally incident light. Both a two- level and a three- 
level A configurations are considered. The band dispersion is characterized as a function 
of the relative position of the atomic resonance and the Bragg frequency corresponding 
to the lattice periodicity. Two main regimes are identified depending on whether the 
two frequency scales are close or well separated. The band dispersion is then used in 
a calculation of the reflection spectra for both semi-infinite and finite systems. These 
results are the starting point to propose a photon energy /z/ter experiment [12], where the 
carrier frequency of a slowly travelling wavepacket can be continuously tuned while fully 
preserving its pulse shape and its coherence. 

The paper is organized as follows. In SeclH we review the Hopfield approach to radiation- 
matter interaction and its application to the case of a lattice of two-level atoms. In Secl3l 
the Hopfield formalism is applied to get predictions for the dispersion of the elementary 
excitations of the system, the so-called polaritons. Reflection spectra for semi-infinite 
and finite ID systems are presented. In SecJH we extend the model to the case of three- 
level atoms in a A level configuration: we implement the dressed atom approach into the 
Hopfield Hamiltonian and we obtain predictions for the polariton dispersion and reflection 
spectra. In SecJH the application of three-level atomic systems as photon energy lifter 
is discussed for realistic parameters of Rb systems. The main experimental issues are 
discussed in SecJHl Conclusions are drawn in SecE 

2 Hopfield approach 

We consider the interaction between light and a collection of two-level atoms trapped in 
the periodic potential of an optical lattice [8,9]. The atoms are assumed to be in a perfect 
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Mott Insulator state where one and only one atom is present at each lattice site [10,11]. 
With a suitable choice of the lattice laser frequency and polarization, the minima of the 
optical potential felt by the different atomic states can be made to coincide and to form 
a cubic lattice of spacing /. Provided trapping at each site is tight enough, the atomic 
center-of-mass motion results frozen in the motional ground state around the potential 
minimum (Lamb-Dicke regime) and no atomic tunneling is allowed between adjacent sites. 
The energies of the ground \g) and excited |e) states of the atom have therefore to include 
the light shift due to the optical potential as well as the zero-point kinetic energy: the 
resulting energy separation is indicated by tiw^. 

Single-atom excitations on the i^^ atom are respectively destroyed and created by the 
operators hi = \g)i (e|j and b] = \e)i{g\i. In a solid-state terminology, these localized exci- 
tations can be seen as an extreme kind of Frenkel excitons where the electronic excitation 
is confined to a single atom or molecule of the crystal. This, in contrast to weakly bound 
Wannier excitons whose electron and hole can be spatially separated by a distance much 
larger than the lattice spacing [4]. 

In addition to the coupling to the transverse e.m. field that we shall discuss at length in 
what follows, excitons are coupled to the longitudinal e.m. field. This Coulomb interaction 
between the dipoles of distinct atoms allows for the transfer of excitation from one site to 
another according to Hamiltonian terms of the form fej bj [14] . To take the most advantage 
of the translational symmetry of the lattice, one can construct creation operators for 
delocalized excitations with a well-defined Bloch wave vector k: 



An integration box of size L is assumed with periodic boundary conditions. = (L/l)^ 
is the total number of atoms in the lattice (which is assumed to fill the whole integration 
box). The Ck, cj^ exciton operators satisfy the following approximate Bosonic commutation 
rules 



where M is the number of excitations present in the system. Excitons therefore behave as 
bosons, at least in the "linear optics" limit M <^ in which the probability of a double 
excitation of the same atom is negligible. Throughout the whole paper, we shall stick to 
this limit. 

Thanks to the translational invariance of the system, the Hamiltonian describing the 
internal atomc dynamics is diagonal in the operators: 



k 

The k dependence of the exciton band ti;o(k) is a consequence of the Coulomb dipole- 
dipole interaction and describes the exciton propagation [21]. As this dependence is quite 
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weak in the present system, we will for simplicity neglect it in the following and take a 
constant value ujq throughout the whole first Brillouin zone (fBz) of the reciprocal lattice. 
In physical terms, the difference between uq and ujq corresponds to the Clausius-Mossotti 
local field correction [14]. 

The quantized transverse electromagnetic field is represented by the e.m. vector potential 
operator [13] 



where ak+g,A and a^+g x are the photon annihilation and creation operators for the different 
modes, labelled by their polarization state £"a=i,2 and their wave vector. This latter sum 
is split into the sum over Bloch wave vectors k spanning over the fBz of the reciprocal 
lattice, and vectors g belonging to the reciprocal lattice. In the cubic lattice geometry 
under consideration here, the reciprocal lattice is itself cubic, with a lattice constant equal 
to 27r// [4]. The free- field Hamiltonian has the usual form 

H field = ^ ^k+gol+^a^+s (5) 

with the vacuum frequency of the photon cuk+g = c |k + g|. c is here the speed of light 
and eo the dielectric constant of vacuum. Throughout the whole paper SI units are in use. 

In addition to the terms describing the internal atomic dynamics Hat and the non- 
interacting field Hfieid, the total Hamiltonian has to include terms Hint that couple the 
matter to the transverse e.m. field. This is most simply done by means of the standard 
minimal coupling replacement [13]. As we are considering optical fields with wavelengths 
much bigger than the atomic size, a dipolar approximation can be performed within Hint- 
Atoms are represented as electric dipoles of dipole moment /i. For the sake of simplicity, 
we restrict to a single polarization state parallel to the dipole moment of the atoms. In 
this way, we can drop the polarization index A and write Hint in the final, compact form: 



Hint = ^ (ck«jc+g + c_kak+g) + h.c. 



+ 



+ Yl ^===(4+g«k+g' + 4+g«lk+g' + h.c). (6) 

k,g,g' V^k+gt^k+g' 

We can easily recognize the different interacting terms: the former is the usual dipole 
exciton-photon coupling, with strength proportional to 



"^ = /^^oa/^, (7) 

and the latter is a "diamagnetic" photon-photon coupling induced by the presence of 
atoms, with a coupling strength proportional to 

2 

m' = — . 8 
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For typical values of the system parameters, the coupling coefficient m is generally small, 
i.e. m/ y/uj^ ^ hujQ, which implies m' /uq <^ m/ 

Thanks to the discrete translational symmetry of the lattice, Bloch momentum is con- 
served by all terms in This means that only states with the same Bloch wave vector 
are mixed by the light-matter interaction. As this set is discrete, radiative decay (i.e. 
spontaneous emission) is forbidden and energy coherently oscillates between the photonic 
and the excitonic degrees of freedom [14]. This observation is crucial in simplifying the 
physical analysis of the system. 

The normal modes of the system, the so-called polaritons are therefore superposition of 
the original photonic and excitonic modes; they are classified by their Bloch momentum 
k, and by a band index n. Thanks to the quadratic structure of the Hamiltonian, the 
polaritonic operators can be obtained by means of a Hopfield-Bogoliubov transformation 
in the general form [14]: 



Provided the correct normalization la^k.np — kk,nP + ^g (l^^k.g.nP — Iz/k.g.kP) 
the polaritonic operators satisfy bosonic commutation rules 



<5k,k' <^n.,n' 



1 is chosen. 



(10) 



The coefficients x,y,w,z are determined by solving the eigenvalue problem associated 
with the commutator 



Q^k,n) 



H 



:iii 



'k,nC^k,n5 

where hQ-^^n is the polariton energy. The solution of this system is equivalent to the 
diagonalization of the following Bogoliubov matrix 



-Mi 




(12) 



with the Hermitian matrix 

/ 



—im 
—im 



hjjj] 



k+g 



2m' 



2m' 



2m' 



v"^k+gt^k+g' 



^k+g' + 



2 m' 
^k+g' 



v 



and the symmetric matrix 



Ml 



(13) 



V^k+g 
— im 



V^k+g 
-2m' 
'^k+g 
-2m' 



-2m' 



V'^k+g'^k+g' 

-2m' 



V^k+g' V'^k+g'^k+g' 



"k+g' 



(14) 
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The first row and column correspond to the matter excitation. For notational simphcity, 
only two photonic modes of wave vectors k + g and k + g' have been explicitely shown 
here, but the matrices are intended to contain rows and columns for each reciprocal lattice 
vector g. 

While the diagonal blocks /Ck of are hermitian, the non-diagonal Ai^ ones break the 
hermiticity of the matrix in the usual sense: the matrix £k is in fact 0-hermitian, in 
the sense that QC^Q = where = diag(l, 1, 1, —1, —1, —1) defines the Bogoliubov 
metric. Physically, this property is related to the fact that the A4 blocks correspond to 
the anti-resonant terms in the Hamiltonian which do not conserve the number of ex- 
citations [15-19]. Thanks to the small value of the light-matter interaction coefficients m 
and m' in the atomic systems under consideration here, most of the physics under inves- 
tigation here can be obtained by neglecting the anti-resonant terms and truncating 
the /Ck matrix to a small number of photonic modes [20,21]. 

3 ID lattice of two-level atoms 
3.1 Band structure 




Figure 1: Polariton dispersion in a ID lattice of two- level atoms. Panel (a): Purely 
excitonic regime, exciton-photon coupling (m/y/uJ^) ^ 8 x 10~'^hu!o, distance between 
resonance and Bragg frequencies ubt—^o ~ Scjq- Panel (b): Mixed exciton-Bragg regime, 
exciton-photon coupling (m/^/cJo) ^ lO^'^hwo, distance between resonance and Bragg 
frequencies ujbt- — ciJo ~ 1.4 x IO^^cuq. The gray regions correspond to the gaps. For the 
sake of clarity, the dipole moment has been exaggerated with respect to actual values of 
atomic systems. 

To get a simple physical understanding of the system, it is useful to concentrate our 
discussion on the simplest case of a ID geometry: most effects related to resonant light- 
matter interaction are in fact indipendent from the dimensionality of the system under 
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consideration [6,20,22]. Interesting discussions of the optical properties of different kinds 
of ID resonant PCs can be found in [7,23-27]. 



Two frequency scales are to be considered: the atomic resonance frequency luq, and the 
Bragg frequency = en /I of the lattice which carries information on the periodicity of 
the lattice and derives from the famous Bragg condition for diffraction spectroscopy in a 
reflection geometry [4]. The width of the frequency region in which states are effectively 
mixed and their splittings are determined by the interacting Hamiltonian ([6]). At leading 
order, this is proportional to mj^/To. Starting from this consideration, two different 
regimes can be distinguished according to the ratio between the difference ^(cuq — ^Br) 
and the characteristic mixing m/ ^uj(Q^Br)- 



The first regime, let's call it purely excitonic regime corresponds to the case when the 
resonance frequency uq and the Bragg frequency UBr are well separated h\u;Br — ^^o\ 3> 
m/^Cij(o,_Br)- An example of polaritonic dispersion for this regime is shown in Fig, 1(a) for 
uq < uobt- For the sake of clarity, a somehow exaggerated value of m has been used in 
the figure. 

On one hand, the exciton and a single photonic mode intersect in the interior of the fBz; 
their mixing results in an anticrossing of the polariton modes with a Rabi splitting equal 
to h^uj = 2m/y/uJ^. Far from this region, the polaritonic modes tend to almost purely 
photon and exciton modes. On the other hand, pair of photon modes intersect at the 
edge of the fBz at frequencies multiple of UBr', their mixing is due to Bragg scattering 
processes on the atomic lattice. 

Gaps open at the edges of the fBz: the lower one (around c^o) is the usual polaritonic 
gap of resonant dielectrics [28], while the upper one (just above ubt) is due to Bragg 
scattering processes. Differently from the bulk case, the former gap extends on both sides 
of uq because of the limited size of the fBz. The latter one is instead located strictly 
above UBr- Its lower edge is exactly at UBr and corresponds to a purely photonic state, 
unaffected by the presence of the atoms which are located at the electric field nodes. As 
usual, the polaritonic density of states vanishes inside the gaps, and radiative propagation 
at these frequencies is forbidden. 

From a formal point of view, the lower gap mainly originates from the exciton-photon 
coupling, while the upper one contains a contribution from the direct photon-photon 
coupling as well. As m? /{fujjBr) ~ ■, all these terms are of the same order. This picture 
is preserved in 2D and 3D, as all other photonic modes that may participate are far away 
in energy. 

The other regime, let's call it mixed exciton- Bragg regime, is characterized by the condition 
fiiujBr — ^o) ^ m/ ^uJ{i^^bt)- Ih this case three modes are simultaneously mixed, namely 
two photon branches (the incoming one at k and the first Bragg diffracted at k ~ 211/1) 
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and the excitonic state. Differently from the previous case, the splittings (again of the 
order of m/ yJuj(Q^Br)) are now located close to the edges of the fBz. As one can see in Fig. 



1 (b) , this results in much wider forbidden gaps of the order of the splitting. In addition 
to this, the "central" band between ojo and ojbt results squeezed and shows a very flat 
dispersion over most of the fBz. 

In the mixing region around the band edge, the upper and lower polaritons are mixture of 
exciton and photons with almost equal weights, while the central polariton band at ujbt 
is mostly photonic. This effect is easily explained in terms of the electric field showing 
nodes at the atomic locations, as in the previous regime. 



3.2 Reflectivity spectra 

The band dispersion introduced in the previous section is a complete description of the 
photon propagation inside the lattice. Most spectroscopic experiments, however, involve 
light beams which are incident onto finite systems and therefore require a description of 
the interfaces between regions of different optical properties, namely the external vacuum 
and the atomic lattice. This allows to calculate crucial properties of the system, such as 
its reflectivity [7,23-26]. 

This can be done by imposing suitable boundary conditions on the electromagnetic fields: 
Maxwell theory imposes in fact the continuity of both the electric field and its spatial 
derivative. In the following, we shall restrict to the most relevant case of a coherent, 
monochromatic excitation at uj. Outside the lattice, we have purely photonic, coherent 
plane waves with wave vectors such that \k\ = uj/c. Inside the lattice, the field propagates 
as coherent polaritonic states of modes such that Qk,n = ^■ 

The electric field is the expectation value of the electric field operator 

E{x) = zXAI J2 (se^"^ - ay^^^) . (15) 
g 

Outside the lattice, the electric field associated with a photonic coherent state of wave 
vector k is given by 

E{x,t) = {coh : ak{t)\E{x) \coh : a,{t)) = m/tt^ e*^'^™*) + h.c. (16) 

Inside the lattice, one is to consider a polaritonic coherent state with Bloch wave vector 
k in the n band satisfying 

dfc,„ \coh : ak,n{t)) = ak,n{t) \coh : ak,n{t))] (17) 
time evolution of such a state under the system Hamiltonian corresponds to 

«fc,n(t) = 6-*"'=-* = e-'^K (18) 
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The electric field of such a state is then: 



E{x,t) = {coh : ak,n(t)\E{x) \coh : 

E 



J{{k+g)x-nk,nt) 



-g.nJ 



+ h.c. 



■ (19) 



The amplitudes at and ak,n are then determined by matching the fields (|T6l) and (fT9l) at 
the interfaces for the given geometry under consideration. 

The plots in Fig. [2] and Fig. [3] show the reflectivity spectra in respectively the purely 
excitonic and the mixed exciton-Bragg regimes. Two geometries will be considered: a 
semi- infinite [panels (a)], and a finite one [panels (b)]. Note that these predictions exactly 
coincide with the ones of semi-classical calculations where matter is described in terms of 
a dielectric polarizability [6,7,22,24-27]. 
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Figure 2: Reflectivity spectra in the purely excitonic regime for (a) a semi-infinite lattice, 
(b) a finite one containing M = 100 cells. Parameters as in Fig. [TJ The gray regions 
correspond to the gaps. 



3.2.1 Semi-infinite geometry 



In the first case. Fig. 2(a) and Fig. 3(a), there is a single interface, dividing the space in 



two semi-infinite regions: vacuum and lattice. 



We first consider the input problem with an incoming and a reflected plane wave in 
the vacuum and a single transmitted polariton Bloch mode in the lattice: the Bloch 
momentum k in the fBz and the band index n of this latter are chosen in order to satisfy 
energy conservation u = Qk,n- The reflection amplitude r^^ is then expressed in terms of 
the electric fleld of the polariton and its spatial derivative as 



,2KWc)(//2) Eil/2)+tic/uj)E'il/2) 
E{l/2)-i{c/uj)E'{l/2y 



(20) 
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Figure 3: Reflectivity spectra in the mixed exciton-Bragg regime for (a) a semi-infinite 
lattice, (b) a finite one containing M = 100 cells. Parameters as in Fig. [H The gray 
regions correspond to the gaps. 



The point x = 1/2 lies on the edge of the elementary cell, i.e. midway between neighbour- 
ing atoms 0. 

The mismatch between the photonic components of the polariton state and the incoming 
wave determines the reflectivity Rin = |?"mP shown in the figures: this is significant around 
the gaps where the incoming wave is strongly mixed with either an exciton [around the 
lower gap in Fig 2(a)| , a photon [around the upper gap in Fig 2(a)| , or both [around both 



gaps of Fig, 3(a)| . In the mixed exciton-Bragg regime, note that the reflectivity remains 



quite large in between the two gaps: the flatter the middle-polariton branch, the higher the 
corresponding reflectivity. In the present semi-inflnite geometry, reflectivity is complete 
inside the gaps where the wave vector becomes imaginary and the fleld inside the lattice 
consists of an evanescent wave. 

The output problem corresponds to two counterpropagating Bloch modes with the same 
energy inside the lattice and a single transmitted plane wave in the external vacuum. In 
this situation, the reflectivity is given by Rout = koutP with 

E4l/2) + ^ic/uJ)E'4l/2) 

E_{l/2) + i{c/uj)EL{l/2y ^ ' 

where the -|- and — signs refer to the propagation versa of the polaritons. 

The system being invariant under time reversal and spatial parity, it is easy to prove 
that the coefficients Wg and i/g entering in the formula (1191) for polariton fleld share the 



^Actual numerical calculations arc performed by truncating the full matrix (|12p to a finite number 
of modes; in order to obtain a smooth convergence, a gaussian cut-off has also been added on higher 
photonic modes. This cut-off physically mimicks finite size atoms. Because of the singularity in the 
field at the atomic position x = 0, evaluating the fields at x = 1/2 rather than at x = ensures faster 
convergence. 
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same phase and can be chosen to be all real. The electric fields of counterpropagating 
polaritons are therefore complex conjugates of each other E_{x) = E^{x)*. Plugging this 
into ([21]) and comparing the result with fl20l) . it is immediate to see that Rin — Rout- 



3.2.2 Finite slab 



Refiectivity spectra are shown in Figs. 2(b) and 3(b) for a finite system containing a 



quite large number M of elementary cells. Note that the present approach is not able to 
calculate the refiectivity inside the gaps where the field consists of an evanescent wave. 
The main difference with respect to the semi- infinite case is the presence of fast oscillations 
on top of the refiectivity spectrum around the main gaps. This can be explained as follows. 

Two interfaces at respectively Xfr = —{{M — 1) + 1/2)/ and Xback = now separate 
three regions of space: the vacuum with the incident and refiected photons, the finite-size 
lattice with counterpropagating polaritons, and again vacuum with now only a transmitted 
photon. The field in the last cell {x G [— //2,//2]) is determined by the output problem 
to be 

Est{x) = E+{x)+routE4x). (22) 

As both E±{x) are Bloch states, the field in the first cell (taking x G [—{{M — 1) + 
1/2)/, -((M - 1) - 1/2)/]) has the simple form 

Est{x) = E+{x + (M - 1)/) e-^'^'^^^-i) + routE^{x + (M - 1)/) e''^'^^'-^\ (23) 

By solving the continuity conditions at the front interface at x = Xfr, we get 

^^^^^ ^ ^2i{u.lc)xr. Est{Xfr) + i{c/uj)E'JXfr) ^^4) 
Est{Xfr) - i{c/uj)E'^i{Xfr) 



Because of the phase factors in 0231) . fast oscillations occur in the refiectivity 0241) due 
to the Fabry- Perot- like interference of Bloch waves which undergo multiple refiections at 
the lattice boundaries. The period Acj of these oscillations is fixed by the group velocity 
"^fc^n ~ dQk n/dk and the total length of the system L = Ml, 



the slower vl^^, the closer the peaks. 



To compare the envelope of this oscillations with the spectrum in the semi-infinite geom- 
etry, we can consider a simplified model where the lattice is replaced by a bulk medium 
of refractive index n. In this case, the refiectivity for a single interface separating vacuum 
and medium is 

fl.„.=(i^)^ (26) 

For a slab of thickness L, the refiectivity is [29] 

= jn-l/n)^ sin^(^nL/c) ^^7) 

4 cos^ {unL /c) + {n + l/nY sin^ {unL / c) 
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Fabry-Perot oscillations are apparent, with a maximum reflectivity at the peaks equal to 

1 2\ 2 

1- n-^^ 



RTa^ = [yt^ ) ■ (28) 

In the limit n — 1, the ratio {BJ^iab /Rmt) 4: this is due to the presence of two 
counterpropagating Bloch modes in the slab as compared to the single propagating mode 
in the semi-inflnite case. This factor 4 provides a good approximation in the lattice case 
as well, as one can easily see in the low-reflectivity tails of the spectra of FigsI2] and El 

The peak at Bragg frequency, which is more related to the photonic components, is 
generally smaller than the resonant one. This difference is dramatically enhanced in the 



purely excitonic regime as you can see in Fig, 2(b) 



4 ID lattice of three-level atoms 

Systems of three-level atoms are of great relevance for both quantum and non-linear 
optics. A most remarkable phenomenon in this respect is the so called Electromagnetically 
Induced Transparency (EIT) [30] which, among other properties, can lead to ultra-slow 
light propagation [31-33] in spectral regions where absorption and reflectivity are also 
very low. This fact makes three-level systems extremely interesting systems to study and 
manipulate light in new regimes. 

We consider the three-level A conflguration shown in FigJU in addition to the previously 
considered two- level scheme, there is a metastable state \m) which is coupled to the excited 
state by a classical (laser) fleld. The |m) l^f) transition is optically forbidden. 




Figure 4: Sketch of the three-level A conflguration. 



The coupling beam dresses the metastable state, so to give a new mixed excitation formed 
by the atom promoted to the metastable state and an extra photon correspondingly 
emitted into the dressing beam: this new state has an energy fujm = f^il^m + ^c), where 
ojm and Uc are the frequencies of respectively the bare metastable state (with respect to 
the ground state \g)) and the dressing laser light [13]. Its detuning from the excited state 
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is then 5c = {uim — ^^o)- The Hamiltonian describing the dressing is 

Hc = ^Yl i^ld. + ^^^"0 ' (29) 

i 

where Qc is the Rabi frequency of the dressing beam and the operators di, d\ destroy and 
create the dressed metastable excitation for the atom i. 

As above, the introduction of the radiation-matter interaction is performed by means of 
the minimal couphng replacement [19]: 

hi-^hi-i^A{k). (30) 

The d operators remain instead unchanged as the metastable state is not directly coupled 
to the quantized field. After the minimal coupling replacement, the dressing Hamiltonian 
has the form: 

Hc = ^Yl {^^^^ + ^^^1 + '^Mh)di - ^^i(/i)4) • (31) 



Note the appearance of terms directly coupling the dressed metastable state to the light 
field. These terms are however of the order of ^/loq (with = mfic/StUo), and 
therefore much smaller than all the other coupling strengths in the system. 

By following the same procedure as above, new excitonic delocalized operators can be 
constructed for the dressed state, and then included in the Bogoliubov matrix: the result- 
ing polaritonic bands and the corresponding reflectivity spectra are shown in Figs. [5](a,b) 
in the purely excitonic regime and for a resonant dressing 5c = 0. Note that only the 
region around the atomic resonance is shown in the figure. 

For the sake of simplicity, we concentrate our attention onto a purely excitonic case cuq ^ 
lobt where the lattice spacing is much smaller than the resonance optical wavelength. The 
signature of the three-level scheme is the presence of a fiat EIT mini-band corresponding 
to the "central" polariton. The most interesting region is around the Raman resonance 
ojk = 'j^m between the quantized radiation and the dressed state. This spectral region can 
be simply explored in terms of the reduced 3x3 Bogoliubov matrix 



^ h6k imj^Juj^ ^ 
—im/y/uJ^ —h6c 
\ / 



(32) 



where only the incident photon (1st row) mode, the excited state (2nd row), and the 
dressed metastable state (3rd row) are considered. Sk = ck — Um is the detuning of the 
photon mode from atomic resonance. As we are considering a small region 6k ^ ujq, we 
have replaced Uk with ujq in the light-matter coupling terms. We have also neglected the 
m' and nic coupling terms as they are much smaller than the other terms. 
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Figure 5: Polariton properties in a ID lattice of three-level atoms. Dispersion (a,d), 
reflectivity (b,e), group velocity (c,f) in the region around resonance. Red, solid line 
is lower polariton; blue, dot-dashed line is central polariton; green, dotted line is upper 
polariton; black, dashed line is photon dispersion in vacuum. Parameters for a / = 100 nm 
lattice of Rb atoms in the purely excitonic regime: ubt — ~ Sc^o? exciton-photon 
coupling m/^/cJo ~ 4 x lO'^^hwo, dressing amplitude 0,c/ujq = 2 x 10^^. Panels (a,b,c): 
u! ^ ujrn = 1^0) = (initial state of the photon lifter). Panels (d,e,f): 6c/ujq = 1.2 x 10~^ 
(final state of the lifter). In the different panels, the blue circles indicate the position of 
the wavepacket to be "lifted" . 
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A noticeable property of ( l32l) exactly at Raman resonance 5^ = is that for any value 
of different parameters there exists an eigenvalue A = corresponding to the central 
polariton. Its eigenvector has a vanishing exciton component and a photonic weight 



so that the group velocity is 



Q2 

= + {Amyh'uo) ' ^^^^ 



" fi2 + (4mV/i2a;o) " ^^^^ 



a significant reduction of the polariton group velocity is then observed as soon as a weak 
dressing amplitude is used HQc <^ '^/v^ (see FigJSt^c)). Furthermore, reflection at the 
lattice interfaces is vanishing in the central region of the central polariton band, which 
allows for easy injection of polaritons in the system (as shown in FigEI^b)). 

A simple model to quantify the width of the region where reflection is low can be developed 
as follows. Defining an effective index for the lattice as Ue/f = ck/Qk = {cko + 6^)/ {cko + 
6kVgr/c), and applying classical reflection result (PHj) . the reflectivity envelope for a semi- 
infinite system turns out to be approximately given by: 

R^nt - —^-T- (35) 

This approximation is accurate at the center of the dip where the reflectivity is small. 
Comparing the two main results (IH^ and of the present section, a trade-off is appar- 
ent: the slower the light, the narrower the reflectivity dip [32]. 



5 Application: photon energy lifter 

Obtaining a coherent and widely tunable frequency conversion of an optical signal is a 
central task in optical telecommunications [34]. Several techniques have been developed 
during the years to perform this operation, but most of them suffer from significant 
limitations in their application range, or are disturbed by spurious effects. 

A recent proposal is based on the so-called dynamic photonic structures (for a review see 
e.g. [35]), whose optical properties are varied in real time while the optical wavepacket 
is propagating inside them. The basic idea of the photon lifter consists in the adiabatic 
shift of the photonic band on which the photon is located. This was orginally proposed 
for solid-state photonic structures [12], but it is interesting to explore the potential of cold 
atom systems to this purpose: the very long coherence times in a Mott insulator state 
and the easy tunability by external electric or magnetic fields makes them very promising 
for this kind of applications. 

As a specific example, we shall consider in what follows a system of ^^Rb atoms, trapped 
in a cubic optical lattice of spacing / = lOOnm. The optical properties are varied by means 
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of an external magnetic field (taken as parallel to the z axis) which shifts the atomic levels 
via the Zeeman effect [36]. 

In the following, we concentrate on the D2 transition at a frequency ujq/2ti = 384 THz [37]. 
As we are interested in substantial shifts, we can concentrate our attention in the high 
field regime {B > 5 x 10'^ G) where the atomic nucleus is decoupled from the electronic 
degrees of freedom, and the energy shift mostly comes from the electronic total angular 
momentum only: AE = ^bQjJzBz-, where /i^ is the Bohr magneton, gj is the Lande 
factor of the considered level and is the z component of the total angular momentum 
of the electron. We use the \ J = 1/2, = =Fl/2) sublevels of the 5^S'i/2 electronic ground 
state as respectively ground \g) and metastable |m) states, and the | J = 3/2, = 1/2) 
sublevel of the 5^P3/2 electronic excited state as excited |e) state. The corresponding 
Lande factors are gj =1/2 = 2 and (77=3/2 = 4/3. The nucleus is not affected by the optical 
process and maintains the same polarization it had in the initial state: in the absolute 
atomic ground state, the nuclear spin is e.g. polarized antiparallel to the electron spin of 
the I J = 1/2, J, = -1/2) state. 

A z polarization is used for the dressing light beam that couples the |m) and the |e) states, 
and a 0"+ one is used to probe the polariton dispersion on the \g) \e) transition. Using 
tabulated values for the electric dipole moment of the D2 transition, the exciton-photon 
coupling ([7]) for the system under consideration is of the order of m/^/cJo ^ 10~^/ia;o. 

To maximize the available time to perform the lifter operation, it is useful to have a very 
slow group velocity, which in turns requires a small dressing amplitude. In the following, 
we shall choose Qc/^^o = 2 x 10^'^. This value Qc/'^'^ ~ 76 MHz corresponds to 10 times 
the radiative linewidth of the D line of Rb atoms. 

The dressing frequency is chosen in a way to have 5c = at the initial value Bin of the 
magnetic field: the corresponding polariton dispersion is the one shown in Figl5](a). The 
light pulse is injected into the system in proximity of the resonant point 6k = 0, where 
the interface reflectivity goes to zero, and injection is most effective (see the circle in 
Fig|5](a)): the width of this dip results from (l35i) to be of the order of 2 x 10~^Ci;o and the 
group velocity is Vgr/c ^ 7 x 10~®, i.e. Vgr ~ 20m/s. 

The magnetic field variation is performed while the light pulse to be shifted is completely 
contained in the lattice and is propagating through an effectively bulk system. As the 
magnetic field is varied in a spatially homogeneous way, the Bloch wave vector is conserved 
during the process. If the field variation is slow enough as compared to the frequency 
difference of neighbouring bands, the polaritons will adiabatically follow the band and 
their frequency at the end of the process will be accordingly shifted (see the circle in 

FigEKd)). 

As an example, we propose to tune the magnetic field from 1 up to 2 T: this results in 
the metastable and excited states being shifted by respectively {6m — 6g)/u}Q = 7.3 x 10~^ 
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and ((5e — 5g)/ujQ = 6.1 x 10^^ with respect to the ground state. For hght initially injected 
in proximity of ujq = the shift of the photon frequency results approximately equal 
to 6m, which amounts to the quite sizeable value 14GHz/T. As the lifter operation is 
based on an adiabatic shift of the polariton dispersion, it completely preserves the pulse 
shape and the coherence properties of the incident wavepacket, both at classical and at 
quantum level. 

6 Experimental issues 

To verify the actual feasibility of such a promising experiment, it is important to mention 
the main practical difficulties that may arise in an actual experiment, and discuss how 
these can be overcome. 

1. We have verified that the transmittivity of the lattice interfaces is close to 1 for 
both the injection and the extraction process (FigsJH^b,e)). The pulse is injected 
into the lattice at a frequency corresponding to the EIT refiectivity dip around 
Raman resonance. The extraction takes place in close proximity of the Raman 
resonance where refiectivity is again very low. This, in spite of the fact we are very 
close to a gap: thanks to the now significant detuning 6c, the metastable state is 
in fact weakly coupled to light, and the corresponding crossing point is displaced 
slightly away from the light line. 

2. In order to have a reasonably long time to vary the magnetic field, we have verified 
that the group velocity of the polariton states involved in the lifter operation is slow. 
Light initially propagates on the EIT slow light branch, which is deformed during 
the lifter operation. At the end, the wavepacket is found on the very fiat region 
below the gap where the group velocity is low. 

3. The wavepacket has to be shorter than the lattice length, still its frequency spectrum 
has to fit in the refiectivity dip at both injection and extraction. A lattice of M = 
1000 cells is able to accomodate pulses with at most Ak > l/(/M) = ksr/i'^M). 
From panel (a), this corresponds to a lower bound on the frequency width of the 
incoming wavepacket Aum = Akvf^ > 5 x 10~^^ loq. One can easily see in panel (b) 
that this frequency spread still fits within the injection window where refiectivity is 
low. The same can be verified on panels (d-e) for the extraction process. 

4. In order for the pulse shape not to be affected, dispersion of the group velocity 
should be small for the wavevector window Ak under examination. Initially, this 
is not a problem, as we are working close to the center of the EIT branch where 
the group velocity has a weak dispersion. The situation can be more critical on 
extraction, because of the strong squeezing of the polariton band in the region just 
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below the gap. The importance of this effect can be reduced by choosing pulses 
initially tuned just above the Raman resonance. 

One major constraint that still exists on the experimental parameters concerns the speed 
at which the magnetic field has to be actually varied. As this has to be done while the 
wavepacket is inside the lattice, a very slow group velocity and a long lattice are required. 
Using values for state-of-the-art Mis. namely L = M I = 100 /im, and Vgr = 20m/s, one 
obtains that one disposes of a time of approximately 5 fis to perform the magnetic field 
variation. This means that a variation of AB = 1 T requires a very large rate of 2 kG/ /is. 

As this can pose serious difficulties in an actual experiment, it is worth briefly exploring 
alternative strategies. An interesting possibility is to further reduce the dressing amplitude 
Qc- As the polariton group velocity is proportional to the square of the dressing amplitude, 
the value fi^ = 2 x 10~^uq used in a recent slow light experiment [33] already leads to a 
group velocity of the order of 20 cm/s which corresponds to an available time of 500 /is. 
In the high-field regime considered here, a photon frequency shift of 1 GHz then requires a 
magnetic field variation of 500 G in 500 fis, a rate routinely used in cold atom experiments. 

It is important to note that the reduction of the dressing amplitude implies a squeezing 
of the reflectivity dip at injection and an enhancement of the dispersion at extraction. 
These, apparently serious problems are overcome thanks to the fact that a reduction in 
the group velocity implies a spatial shortening of the pulse in the lattice, and therefore a 
reduced frequency spread for a given length. 

7 Conclusions 

In conclusion, we have developed a fully quantum description of radiation-matter interac- 
tion in a gas of ultracold atoms trapped in the Mott insulator phase of an optical lattice. 
The coherent interaction between photons and atomic excitations gives rise to new, mixed 
polaritonic excitations. 

In the case of two-level atoms, two different regimes are identified. In the purely excitonic 
regime, where the atomic resonance is far from the Bragg frequency of the lattice, two 
gaps appear in the energy spectrum. In the mixed exciton-Bragg regime, the interplay of 
the atomic resonance and the lattice periodicity enhances the gap amplitudes and gives 
rise to a fiat band between them. The consequences of the polariton dispersion on the 
reflection properties of finite lattices have been investigated. 

The theory is then extended to a system of three-level atoms in a A configuration. The 
slow-light band which joins weak reflection to ultra-slow group velocity turns out to be 
the most promising in view of applications as a photon lifter, i.e. a device able to shift the 
carrier frequency of an optical pulse without affecting its shape nor spoiling its coherence 
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properties. Advantages and disadvantages of using atomic g photon lifters are 

discussed. 

Future work will address the application of slow as well as stopped polaritons as non- 
destructive probes of the microscopic properties of ultracold gases, e.g. the behaviour of 
matter wave coherence across the superfiuid to Mott insulator transition. 
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